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Abstract 



In this work we consider the oscillations and associated resonance of a 
2D Bose-Einstein condensate under periodic and random modulations of the 
atomic scattering length. For random oscillations of the trap potential and of 
the atomic scattering length we are able to calculate the mean growth rate for 
the width of the condensate. The results obtained from the reduced ODE's 
for oscillations of the width of condensate are compared with the numerical 
simulations of the full 2D Gross-Pitaevskii equation with modulated in time 
coefficients. 
PACS: 03.75.Fi;05.30.Jp 

Key words: Bose-Einstein condensate, Gross-Pitaevski equation, resonance, random pertur- 
bation, trap potential, scattering length 



* Corresponding author. Fax: -h 7-712-35-42-91 
E-mail address: fatkh@physic.uzsci.net (F.Kh. Abdullaev) 

1 



I. INTRODUCTION 



The problem of the oscillations of two and three dimensional trapped Bose-Einstein con- 
densates (BEC) have attracted a great deal of recent attention [0], particularly the problem 
where the trapping potential is temporally modulated. The main physical motivation for 
considering such a problem is the experiments of Jin, et. al. |0, which observed resonant 
response of the BEC oscillations to periodically perturbed traps. Some theoretical investiga- 
tions of the problem have included the work of Castin and Dum ||^, in which the resonances 
in the width oscillations have been analyzed using scaling theory; the work of Garcia in 
which it was argued that the oscillations of 2D condensate can be described by the para- 
metric resonances of linear oscillator; and the papers by Pitaevski and Rosch |15|,[16|1 and 
Kagan et. al. [^], where it was shown that a 2D condensate exhibits a harmonic mode with 
frequency 2ujq when the harmonic trap is driven with frequency ujq. Note that the related 
problem of a soliton interacting with an impurity in 2D molecular crystals had been studied 
earlier by Gadidei et al. 0. 

The excitation of resonant oscillations is also interesting as a possible mechanism for 
the stabilization of a BEC with attractive interactions. Such a possibility has been demon- 
strated for the BEC with negative atomic scattering length under temporally modulated 
trap potential if the surface mode is exited ||^. In this work a variational approach is used, 
and the trap was anisotropic. 

It is also interesting to consider the response of a condensate to temporally periodic 
variations in the atomic scattering length. Such a variation of the atomic scattering length 
can be achieved experimentally by varying the magnetic field [|ri|] or using optically induced 
Feshbach resonances 0]. The case where the scattering length varies monotonically in 
time, in particular the case where the sign changes from repulsive to attractive, has been 
studied recently by Dalfovo et. al and Fedichev, et. al. 0. Similarly the effect of a 
periodic variation of atomic scattering length on the tunneling between two condensates in 



the double-well trap, a resonant tunneling, has been studied by AbduUaev and Krankel 10 
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These observations suggest the possibihty of new phenomenon in Bose-Einstein conden- 
sate width oscillations when the atomic scattering length is allowed to fluctuate about 
some mean value. In this paper we study the influence of both periodic and random fluc- 
tuations of the atomic scattering length on the condensate dynamics. As noted in [6], at 
temperatures T ^ ng, g = A'Kh^as/m^ (where n is the gas density, is the atomic scattering 
length, m is the atom mass), only the condensate evolution is pronounced. The perturbation 
of the thermal cloud is small at these temperatures, and in this case the eigenfrequencies of 
small oscillations of the thermal cloud are close to those of the condensate. 

The specific problem that we consider is the 2D Gross-Pitaevskii equation with periodic 
or random variations in the coefficient of the nonlinear term. This problem is interesting 
not only for BEG but also for nonlinear optics too - for example the dynamics of optical 



beams in nonlinear layered waveguides ||T3|. 

The structure of the remainder of the article is follows: In Sect. 2 we describe the model 
and derive the reduced ODE system to describe the dynamics of the condensate width. 
In Sect. 3 we analyze the effect of time-period perturbations on the width dynamics using 
the action-angle variables for the reduced ODEs. We explore the resonances in condensate 
oscillations in Sect. 4, and the dynamics of width oscillations under fluctuating trap potential 
and atomic scattering length in last section. 



II. THE REDUCED ODE MODEL 

It is well-known that wavefunction for a 2D Bose-Einstein condensate in a trap potential 
V{r) is described by the Gross-Pitaevskii equation: 

iHt = + Vtr{r)^ + g{t)\iPW. (1) 

Here Vn^r) = muP'r^ is the trap potential and g(f) = AnTi^as/m, with being the atomic 
scattering length. We will assume that the time-dependent scattering length is constant to 
leading order, with 

as = ao(l + eoW)- 
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In this paper we analyze the case where eo(t), representing the small fluctuations, is a 
periodic function of time as well as the case where eo(t) is a mean zero white noise random 
process: < eo(t)eo(t') >= 2a'^6{t - t'). 



Under the scaling r = tuj,x' = xyh/{2muj), u = ip^^jATrao/uj the Gross-Pitaevskii equa- 
tion can be written in the dimensionless form 

iUr + Am u + s\u\'^u = 0, (2) 

where s > 0, s < correspond to attractive or repulsive interactions between atoms in the 
condensate respectively. 

A number of methods can be employed to investigate the dynamics of the Bose-Einstein 
condensate under temporal variations of the atomic scattering length. One of the simplest to 
apply is the averaged Lagrangian approach. According to this method we take the Gaussian 
anzatz for the fleld [[r^,|l5|fl . 



uir, t) = Air) exp(-^^ + '-^^ + z0(r)), (3) 



To derive the equations for the wavepacket parameters y4(r), a(r), 6(r), 0(r) one should 
calculate the averaged Lagrangian 



L(r) = J rdrL(r,T). 



For the Gaussian ansatz in Eq.(|^) the averaged Lagrangian is given by 



L = -lA'a^iaX + 20. + ^ + 2a262 + - ^ill^). (4) 
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The Euler-Lagrange equations for the functional L{t) lead to the following equations for 
the the width a(r) and chirp &(r). 



Qr = 2ab 

2 - s(t)N 

br = ^-26^-1, (5) 



with N=^J IwprfV = It is easy to see that the above is a Hamiltonian system, with 
the width a playing the role of a position variable, and the chirp b the conjugate momentum. 



One can eliminate b from the above system which leads to the following evolution equation 
for a: 

2 - sit)N 

au + a = , (6) 

where t = \P1t. For the attractive problem with constant atomic scattering length s{t) = 

1 the variational approach predicts the critical threshold for collapse Nc = 2, compared with 
the the exact value Nc = 1.862. 

We comment that, if the initial condition is close to that of the ground state (Townes) 
soliton, \N — Nc\ ^ e 1, then the modulation theory of Fibich and Papanicolaou [0] leads 
naturally to a modulation equation of the same form as the above, where again the constant 

2 is replaced by the critical value 1.862. The qualitative types of behavior exhibited by the 
above equation do not, of course, depend on the exact values of the constants, since all such 
constants can be scaled out. 

It is useful to rewrite Eq.(|^) in the form 

au + a= — {7) 

where Q{t) = Q{1 + e(t)), Q = 2 — N,e = Neq/ (2 — N). In the next section we construct 
the action-angle variables for the above problem. 

III. PERTURBATION THEORY IN THE ACTION- ANGLE VARIABLES 

We begin our analysis of the Eq.(|^) by constructing the action-angle variables. Of course 
any two dimensional Hamiltonian system can always be reduced to quadrature, but the above 
system is particularly nice because the action-angle variables can be expressed in terms of 
elementary functions. First we note that Eq.(^ is Hamiltonian: 

i/(a„a) = | + f/(a),f/(a) = ^ + ^. (8) 

For < 2, Q > and the motion of the effective particle is bounded, and bounded 
away from due to the "angular momentum" barrier If the energy is E then the width 
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oscillates between Omm = y E — \JE'^ — Q and amax = y E + \JE'^ — Q. The minimum of 
the potential \J occurs at = Q^^'^-, with a minimum energy oi E = V^. = \fQ- When 

> 2, Q < there is no local minimum, and there exist solutions for which a — 0, 
corresponding to collapse of a condensate (see Fig.|I]). 

Since the Hamiltonian is conserved, H = E, we have 



da 



dt 



^2E - - Qa-'^ 
which can be integrated up to the following solution for Q > : 



a(t) = yE + ^E2-gsin(2t + V^o)- (9) 
The action variable for a Hamiltonian system is defined to be 

J = ^fpdq, (10) 

where q is the position variable and p the conjugate momentum. For the above oscillator 



the position variable is a, with momentum p = ^y2Ea'^ — ~ Q/a. The integration is easily 
done via contour integration in the complex a plane leading to the following expression for 
the action vaiable: 

J=l{E-^). (11) 

The bottom of the potential well, E = a/Q, corresponds to J = 0. Since the energy is linear 
in the action, E = 2J + \fQ, the frequency of the unperturbed oscillations is constant 

This can, of course, also be seen directly from the solution given in Eq. (^. 
The Hamiltonian for the perturbed problem is given by 

E = pV2 + U{a) + e{t)V = p" /2 + U{a) + ^ 
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, with V{a) the perturbation Hamiltonian. In the action-angle coordinates J, 9 the pertur- 
bation becomes 

Qy/E''-Qcos{e) 
and the perturbed equations of motion are given by 



IV. RESONANCES OF BEC OSCILLATIONS UNDER PERIODICALLY 
VARYING ATOMIC SCATTERING LENGTH 

In this section we consider perturbations e{t) which are periodic in time. For simplicity 
we discuss the case e{t) — sin(Jli), though this can easily be generalized. To analyze the 
resonances of a BEC under such perturbations we use the multiscale expansion method. We 
introduce the slow time T = et and assume a multiple-scales ansatz of the form 

e = e^^\t,T) + ee^^\t,T) + (15) 

J = J(°) {t, T) + ej(^) + ej(^) {t, T) + ... (16) 

The solution to the evolution on the fast scale is obviously 9^^^ — 2t + J — constant. At 
the next order we find the following evolutions for the action J and the slow angle $ on the 
slow scale: 

where H, the slow Hamiltonian, is given by 

sin{Qt)dt 



277 Jo (2J + QV2 + 2 JV2( J + gi/2)i/2 sin(2i + $)) ' 
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It is clear that we have a resonant response whenever uj = 2n, and the period of the 
perturbation is commensurate with the natural period of the condensate. Note that, even 
though this is a nonlinear oscillator the period is independent of the amplitude. The effect of 
this is that there is no "detuning" from the resonant frequency as the amplitude increases. 
Of course when J becomes large and the width oscillations become significant we no longer 
expect the averaged Lagrangian ODE's to provide an accurate description of the dynamics 
of the condensate width. 

The first (nontrivial) resonance occurs for f2 = 2. In this case the slow Hamiltonian is 
given by 

^' ' 2(J + gV2)i/2gi/2- ^^'') 

and the equations of motion are 

dJ dH sin($)Ji/2 



dt 9$ 2(J + Qi/2)i/2gi/2' 



(20) 
(21) 



dt dJ 4JV2(j + gi/2)3/2- 

The phase plane for the oscillator is depicted in Fig. (^). 

Note that the line $ = — 7r/2 is invariant under the dynamics, and along this line J 
evolves according to 

dJ 

IE " 2(J + QV2)i/2gi/2- (22) 

This solution corresponds to a resonant driving, where the variations in the scattering length 
reinforce the width oscillations. It is clear that for large t the action grows linearly, J cx 
t/2Q^/'^ and thus amplitude of oscillations grows like a ~ \/t. It is also clear from the 
phase portrait that all orbits for which $ 7^ 7r/2 are asymptotic to the invariant manifold 
$ = — 7r/2, so for generic initial conditions one expects that the width will grow like a ~ y/t. 
There is, of course, also the solution $ = 7r/2 in which the variations in the scattering length 
are anti-resonant with the variations in the width of the condensate, and act to damp the 
width oscillations. It is easy to see that when $ = 7r/2 the action J goes to zero in finite 
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time. Since generic initial conditions are asymptotic to $ = —ti/2 the $ = 7r/2 solutions 
are unlikely to be observed experimentally, though it is possible that they could be realized 
with some kind of control. 

Next resonance occurs at = 4. In this case the Hamiltonian is given by 

_ Jsin(2$) 

''"(J + gV2)gi/2- ^^^^ 

Again it is easy to see that the line $ = 7r/2 is invariant under the dynamics, corresponding 
to width oscillations which are in phase with the variation of the scattering length. Since 
the action evolves according to 

dl ^ 2Jcos(2$) 

dT {J + Q^I^)Qy^ ^ ^ 

when $ = 7r/2 we have that the growth of the action J is asymptotically linear - J ^ 2t/Q^/'^. 

A. Numerical simulations 

We have conducted some numerical simulations to test the validity of the calculations 
of the last section. We discretize the problem in the standard way, with time step At and 
spatial step /i, so the approximates u{jh, kAt). More specifically we approximate Eq.(H) 
with the following second order accuracte semi-implicit Crank-Nicholson scheme. 



At 



■^HXl - ^'-D + «i - ^'-i)] + Ii'y " + <t))s\u';\']{u'; + u]^'), (25) 

where e(t) is the perturbation term. In the numerical simulations, as in the analysis, the 
perturbation was chosen to be e{t) = esm{Qt). Eq. (|25| ) represents a tridiagonal set of 
equations for unknowns u^-l, u^^'^ and u^^l [j = 1, 2...(A^ — 1)] in a lattice of points, with 
the values of Uq~^^ and being determined from the boundary conditions |^ |r=o= and 
u{r) |r-=oo— > 0. 
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The set of algebraic equations (^) is solved by the vectorial sweep method. In actual 
calculations the typical space step h ranged from 0.01 to 0.005 and time step At from 0.005 
to 0.002 depending on the closeness of Q to the points of resonance. 

The first experiment, shown in Fig. depicts the solution of the Gross-Pitaevski equa- 
tion when the driving frequency is slightly off resonance: the natural frequency of the width 
oscillations is Qq = 2, while the perturbation has frequency Q = 1.9 and amplitude e = 0.1. 
In this plot, as well as all subsequent plots, the solid line represents the solution to the 
full PDE while the dotted line represents the solution to the reduced ODEs. We observe 
oscillations with beats, representing the superposition of low and high frequency oscillations. 
From simulations we seen the good agreement between the full PDE and the reduced ODE 
model. 

Fig.| represents numerical simulations of the ODE and PDE models at the resonance 
point Q = Qq = 2. The agreement is quite good for the period of the width oscillations, 
although there is clearly some discrepancy in the actual value of oscillation amplitude. Note 
that the same phenomenon has been observed in numerical simulations of the resonances 
in condensate oscillations under periodically varying trap potential The graph in Fig.|^ 
depicts the energy versus time for the ODE (dotted line) and PDE (solid line) simulations 
for the same values of parameters as in Fig.^. The agreement between simulations of the full 
PDE and ODE is very good for time less than 30 or so, though for times between 30 and 40 
the oscillations in the energy of the PDE are smaller than the analogous oscillations of the 
ODE, probably due to radiative damping. As was argued earlier the variational approach 
is unlikely to be valid when the oscillations have large amplitude and other effects, such as 
radiative damping, become important. 

Fig.^ depicts the oscillations of the square width of the condensate at the 2 : 1 resonance 
Q = 2Qq, where the frequency of the perturbation is twice the natural period of the width 
oscillation. As in previous case the frequencies agree very well, but the amplitude of oscil- 
lations is larger for PDE in comparison with ODE. This discrepancy grows with time. In 
this case it seems clear that the Gaussian anzatz does not correctly capture the behavior of 
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the underlying PDE. The rate of growth of the energy, which is not depicted, is still linear, 
in agreement with the theoretical estimates. 



V. EVOLUTION OF BEC UNDER RANDOM FLUCTUATIONS 

The action-angle formulation also provides a suitable framework for the analysis of con- 
densate oscillations under random (white-noise) perturbations. One interesting quantity 
which can be calculated is the mean time to achieve a given distortion under a stochastic 
perturbation. From the ODE point of view this problem is equivalent to the problem of the 
mean time to achieve the given level of the oscillations amplitude in the effective potential. 

The starting point for this calculation is Eqn(|l3|,|l^), where e{t) is taken to be a white 
noise process e{t)dt = dB, with dB the increment of a Brownian motion B. This pair 
of stochastic ordinary differential equations has an associated Fokker-Planck equation for 
P( J, 6, t), the probability of having values J, 9 at time t. In the weak noise limit one can do 
a straightforward multiple scale expansion on this Fokker-Planck equation. Upon doing so 
and averaging over the fast angle variable one is lead to the following equation for P{J,t): 

BP r) BP 

— = a^ — (A(J)— ) P(J,0) = <5(J-Jo), (26) 
where A[J) the average diffusivity for diffusion across energy levels is given by 

1 r2-K 

A{J) = - Ve'dO. 
Zn Jo 



For details of this calculation see the paper of AbduUaev, et. al [|T8[ This one-dimensional 
diffusion can be analyzed in some detail. For instance the mean time to reach action J 
staring from Jq is given by 

Note that if A{J) grows sufficiently rapidly, so that the above integral converges, then the 
mean time to random walk to J = oo is actually finite. 
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A. Fluctuating trap potential 



One physically interesting random perturbation is when the strength of the trap is allowed 
to fluctuate randomly. Since the trapping potential is imposed optically, by a laser, it is 
important to take into account the fluctuations of the effective trap potential due to the 
fluctuations of the laser fleld intensity. The fluctuations of the laser fleld intensity lead to 
random variations of the frequency of the effective harmonic trap u"^ = uJq{1 + e(t)), where 
e{t) is the white noise process |jl9|. The function e{t) should be thought of as the fluctuations 
of the laser intensity around its mean value Eq: 

E{t) - Eo 



e{t) 



En 



This problem was studied in pO[ using a moment expansion, though without numerical 
simulations of the full stochastic GP equation. Here we consider this problem both as an 
illustration of effectiveness of our technique, and to provide numerics for this problem. In 
this case the perturbation Hamiltonian is 



V{a) = -{E + ^E^-Q sine). 
Then from Eqn.(0) the diffusivity is given by 

A(J) = "'^ + v^. 

Calculating the mean exit time < t > we flnd the energy as a function of the mean exit 
time 

E = 2^{e^'<'>/'-^-). (28) 

While the mean exit time is not the same as the physical time in the limit of weak noise 
we expect that, to leading order in the noise parameter, the mean energy should have the 
same dependence on t. In Fig.|^ we present the results of comparison of the theory and 
the numerical simulations of the full stochastic NLS equation (0) for a = 0.04, N = 1.0. 
For the PDE we have averaged over 50 realizations, and for ODE we have averaged over 
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1000 realizations. The figure shows good agreement between predictions of the theory and 
numerical simulations. For large times the solution to the ODE overestimates the actual 
energy. This discrepancy is probably due to effects such as dispersive radiation, excitation 
of higher modes, etc. which are difficult to incorporate into the variational ansatz. 

From the above we can easily estimate the mean exit time from the potential. In optical 
traps the depth of the effective potential is Uq = 2Er = 2/i^A;^/2m, with k = 27r/A, where 
A is the wavelength. The condensate lifetime can be estimated as the time when become 
E Uq. Then 

Substituting the typical values for the laser intensity fluctuations we find that the mean 
exit time is of the order of seconds, tBEC ^ sec. 

B. Fluctuating atomic scattering length 

In the case where the atomic scattering length is allowed to fluctuate we find that the 
perturbation Hamiltonian is given by 

Q^E'-Qcos{e) 
^ 2{E + ^E' -Qsm{e)f' 

From this is follows that the the effective diffusivity is 

and expressing the energy E in terms of the action J we find 

J(2J>vWWQ) 
^^^^ " ^ • 

Note that the expected time to random walk to infinity is actually finite, since A{J) oc J^. 

Substituting this expression for into Eq. (l27|) we obtain for the time to pass from the bottom 

of the potential well where the action Jq = to the state with the action J 

j^Vme^ , = ^. (30) 

2-ey ^ 2 ' ^ ^ 
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or in the terms of the total energy 

E=^. (31) 

2-ey ^ ' 

From the above it is easy to see that the expected time for the width of the condensate to 
grow to infinity - the mean time for the condensate to break up - is given by 

* 21n2 , , 

t* = 32 

In Fig. ||we compare the theoretical expression Eq.(0) with numerical simulations of the 
full stochastic GP equation with fiuctuating scattering length in the case where a = 0.04 
and = 1.0. The energy has been numerically calculated for the localized part of the 
condensate wavefunction. We observe good agreement between the theory and numerical 
simulations for t < 80. 



VI. CONCLUSIONS 

We have considered the oscillations of a 2D BEG with radial symmetry under periodic 
and random modulations of the atomic scattering length, as well as random fiuctuations of 
trap potential. We have calculated the position of resonances and the energy growth using 
a reduced ODE for the condensate width. We have also confirmed the analytical predictions 
with numerical simulations of the 2D Gross-Pitaevskii equation. In the resonant case the 
frequency of oscillations agrees very well with the predictions of ODE, though the amplitude 
shows some discrepancy at large times that is probably due to the approximate nature of 
the variational approach. For the random modulations of the trap potential and the atomic 
scattering length we have calculated the mean exit time corresponding to the time for the 
amplitude of oscillations to exceed a given value and estimated the magnitude of the escape 
time for real experiments with 2D Bose-Einstein condensates. 
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FIGURES 



FIG. 1. Effective potential to describe evolution of the BEG width. 

FIG. 2. Phase plane for the pertubation Hamiltonian at the resonance Q = 2. 

FIG. 3. Oscillations of the condensate width for AO,/Q = 0.1. The solid line is the variational 
approximation, the dotted line the numerical simulations of 2D Gross-Pitaevskii equation (^. 

FIG. 4. The resonant oscillations when = 2. 

FIG. 5. The growth of the energy in the resonant point. 

FIG. 6. The oscillations of width at the second resonance = 4 

FIG. 7. Fluctuating trap: the growth of the energy when a = 0.04 and N = 1.0. Solid, dot 
and dash lines are for the solution of PDF, ODE and Fokker-Planck equation respectively 

FIG. 8. Fluctuating scattering length: the growth of the energy when a = 0.04 and N = 1.0. 
Solid, dot and dash lines represent the solutions of PDE, ODE and Fokker-Planck equation respec- 
tively 
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